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Coarsening systems under uniform shear display a long time regime characterized by the presence 
of highly stretched and thin domains. The question then arises whether thermal fluctuations may 
actually destroy this layered structure. To address this problem in the case of non-conserved dy- 
namics we study an anisotropic version of the Burgers equation, constructed to describe thermal 
fluctuations of an interface in the presence of a uniform shear flow. As a result, we find that stretched 
domains are only marginally stable against thermal fluctuations in d = 2, whereas they are stable 
in d = 3. 



The dynamics of phase separation in systems quenched 
below their critical temperature is substantially modified 
when an external shear is applied. This is true both for 
the case of conserved dynamics (spinodal decomposition 
in binary fluids), and non-conserved dynamics (coarsen- 
ing in Ising spin systems). In the simplest theoretical case 
where the sample is confined between two boundaries 
moving at a constant relative velocity, intuition suggests 
that the growing domain structure becomes anisotropic, 
with domains highly stretched in the direction of the flow. 
This fact is now well established by many experimental 
101 and theoretical investigations [^|-^ : it has been shown 
that the growth of the domains is enhanced along the 
flow direction, while it is less clear whether the trans- 
verse growth rate is unaffected by the shear or depressed 
by it. In either case, the long time effect of the shear 
is to form a structure of very long and (relatively) thin 
domains. 

A natural question in this context is to what extent this 
layered structure is stable against thermal fluctuations: 
long and thin domains may develop fluctuations, trans- 
verse to their main axis, that grow large enough to break 
them. Similar stretching-and-breaking mechanisms have 
been proposed before , in particular to argue that the 
anisotropic domain growth may eventually reach a steady 
state at late times However, the effect of the shear is 
not only to stretch the domains, but also to smooth their 
surfaces, and the net result of the two effects is difficult 
to predict. In order to clarify this problem a suitable 
model for the growth of a surface under shear must be 
analyzed. 

Unfortunately, not many analytic studies of domain 
growth under shear exist. In particular, conserved dy- 
namics and spinodal decomposition have been analyti- 
cally considered mainly in the limit of infinite dimension 
of the order parameter |3|, where no domain interfaces 
are present. The situation is better in the simpler case of 
non-conserved dynamics, where zero-temperature coars- 
ening under shear has recently been studied analytically 
Q . Our task is therefore to study the effect of thermal 
fluctuations in the case of non-conserved dynamics under 
shear. To this end we will introduce a stochastic equation 
for a scalar field h{x, t) representing the height of a fiuctu- 
ating interface above its flat ground state. This equation 



will take into account the smoothing effect of the shear 
flow, which in fact gives origin to a non-linearity of the 
Burgers type. 

Consider an interface separating two regions with op- 
posite order parameter. When the temperature is small 
but non-zero, there will be some fluctuations, resulting in 
a non-flat proflle parametrized by the height h{x,t). In 
the following, we will indicate with d the total dimension 
of the space, and with d' the dimension of the substrate 
spanned by x, that is d = d' + 1. In the absence of 
shear, the free energy cost of a non-flat profile is given 
by = {a/2) J dx (V/i)^, where a is the surface tension, 
and the corresponding Langevin equation coincides with 
the standard Edwards- Wilkinson (EW) growth equation 

i, 

dth^iyW^h + S, , (1) 

where v is a diffusion coefficient (equal to a divided by a 
kinetic coefficient) and £,{x,t) is a delta-correlated noise, 

{^{x, t)i{x',t')) =DS{x~ x')S{t - t') . (2) 

The noise strength D is proportional to the temperature. 
As is well known, due to its linearity a simple scaling 
analysis of the EW equation gives the exact critical ex- 
ponents. Under the rescaling x — > bx, t — > b^°t and 
h b^°h, we have i/ b''°~'^v and D 5^o-2xo-<i+i£) 
By imposing scale invariance, we obtain zq = 2 and 
Xo = (3 — d)/2. Even though the EW equation is not 
suitable for describing interface ffuctuations when a shear 
is present, it is interesting to see what the EW exponents 
would predict for the stability of the stretched domains. 
For d = 2 the investigation of [Q gives two length scales, 
L||(t) = 0{t) and L± = 0{1) (up to logarithmic factors). 
Within the context of the EW equation the transverse 
ffuctuations will grow as 

hr^t^'>/'''F{t/L^°), (3) 

where is a scaling function with the limiting forms 
F(0) — const, F{s) ^ s^xo/zo fQj. g Jj^ writing 

down (||) we have assumed that L\\ is a fixed length scale, 
while actually it is growing with time. The interpretation 
of (||) is, however, simple. If i ^ the domains are 
coarsening faster than the interfacial fluctuations, and we 



1 



can effectively set Ly to infinity. Then h ~ t^o/^o ^ Qn the 
other hand, if t ^ L^°, the fluctuations are coarsening 
faster than the domains and eventually equilibrate on the 
scale Lii. This corresponds to the large argument limit 



'II . 

of the scaling variable in 
these two limits gives 



giving h L y . Combining 



(4) 



In d = 2, this gives h - i^/^, while Lj_ = 0(1). Thus, 
thermal fluctuations would eventually disrupt the do- 
main structure and destroy the coarsening state. In 
d = 3, on the other hand, it has been found that 
L||(i) = 0(t3/2) j^j^j 2._L(t) 0{t^/^). In this case, 
Xo — (logs), so conventional EW thermal roughening 
would not destroy the domains, since h <ti L±. As we 
shall see, the inclusion of shear in equation (|l|) will mod- 
ify these results. In particular, the instability of the do- 
mains in d = 2 will be reduced, while stability in d = 3 
will be confirmed. 

When a shear flow is present the EW equation must 
be modified. We will consider a standard shear velocity 
profile u, with flow along the x direction and shear gra- 
dient perpendicular to the surface. Let us label this last 
direction by z, such that. 



(5) 



where 7 is the shear rate. In this way we break the sym- 
metry between the x direction and the remaining (d' — 1) 
directions x± . Note that the growth is orthogonal to the 
shear flow: this is what happens to the domain walls in 
the long-time limit of a coarsening process under shear. 
The field h{x,t) is now dragged in the x direction by an 
amount proportional to h itself, that is dt ^ dt + jhdx- 
The correct equation for h(x, t) therefore becomes. 



dfh + ^hdxh = i^xdxxh + i'±\7\h + ^ 



(6) 



where we have introduced separate diffusion constants, 
i/x and i/±. For d' = 1 (d = 2) equation (|^) is nothing 
other than the Burgers equation |^ . This can be mapped 
onto the KPZ equation |^ with space-correlated noise, 
via the transformation h — dxh, i.e. hdh ^ ^d{dh)^, 
yielding the KPZ equation for h. Such a case has been 
studied in pO| |. In generic dimension equation (|^) is 
anisotropic and was first introduced in ||Tl[| , in the context 
of a model of sandpiles, and further studied in A 
standard method for the analysis of stochastic non-linear 
equations is the dynamic rcnormalization group (RG) ap- 
proach, first used in this context in Here, we will 
briefly review those RG results for equation which 
are most relevant for our purpose. We start the analysis 
of equation (H) by finding the bare scaling dimensions of 
the parameters. Under the anisotropic rescaling, 

X ^ bx , — » b'^xi^ , t ¥'t , h h^h , (7) 

we obtain, 



Ux t>" ~ Vx b"" ,7 



b^ 



(8) 



Inserting the EW exponents xo and zq found above 
for 7 = into the scaling equation for 7 gives: 7 — >• 
fj{5-d)/2 showing that dc = 5 is the critical dimension 
below which the nonlinearity becomes relevant. Thus, 
for d < 5 the EW exponents are no longer correct and an 
RG approach becomes necessary. 

Burgers equation is known to be invariant under a 
Galilean transformation. In the present context, we have 
an equivalent symmetry, namely the invariance of equa- 
tion (^ under a coordinate transformation which pre- 
serves the original form of the shear flow. Indeed, from 
(^ we have x ~ 'yht, and therefore if we vertically shift 
the interface, h ^ h + ho, the equation is invariant pro- 
vided that we make the transformation x — > a; + 7/19^. 
Basically, this is just translational invariance in the z di- 
rection. This symmetry is important, because it implies 
that 7 cannot be perturbatively corrected in a RG anal- 
ysis and therefore we must set to zero its bare scaling 
dimension. In this way a relation between dynamic and 
growth exponents is obtained, x + ^ = 1- The remark- 
able feature of equation (H) is that there are two further 
parameters whose bare scaling dimensions are not per- 
turbatively changed, namely i>± and D. This happens 
because the 7-vertex in equation (^) carries in Fourier 
space a factor kx that can generate neither terms 0{k'^) 
in the self- energy, contributing to a rcnormalization of 
nor terms 0(1) in the renormalized noise, contribut- 
ing to the rcnormalization of D |l^ . We can therefore set 
to zero the bare scaling dimensions of i^^ and D, obtain- 
ing in this way two extra relations among the fixed point 
exponents: z — and z = 2% + 1 + (d — 2)^. Therefore, 
for d < 5, we have flit]. 
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From (^ and (|^) we can infer the scaling form for 
the correlation function in Fourier space, defined by 
{h{k, uj)h{k', Lu')) = C{k, Lu)S{k + k')S{uj + Lu'). We have. 



(10) 



where / is a scaling function. 

It is useful to rescale coordinates and field in order to 
identify the effective coupling constant of equation (|^). 
The change of variables, 

h (d/uI-'^/^u''/^-') h, x^ ^ (^i/^.)^/' xj., 

i^(iK)t, ^^[D4/yu'/'-'y\, (11) 

amounts to setting i>x = — D — I and replacing the 
vertex 7 by the effective vertex 
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2-(i/4 d/i'-l/2 



(12) 



As expected, for d — dc — 5 this quantity is dimen- 
sionless. Using standard RG methods (see, for example, 
p3| , pl]JlO| ) it is possible to obtain the one-loop RG flow 
equation for the effective coupling constant U = C7^, 
where c = T (3 - d/2) /[8(47r)'*/2-i]^ j,^^ equation 
reads: 



dU 

~di 



(13) 



The trivial fixed point, U = 0, becomes stable above 
the critical dimension dc = 5, giving the EW exponents 
zq, xoj corresponding to ^ = Co = 1- No other physical 
{U > 0) fixed point exists in this phase. The effect of 
the shear on thermal roughening is therefore negligible 
above dimension dc — 5 and domain wall fluctuations 
are isotropic (in the d'-dimensional subspace parallel to 
the mean orientation of the wall). On the other hand, 
in dimension d < 5 there is an attractive flxed point at 
[/ = 2e/9 + O(e^), with e = d^ - d. The effective cou- 
pling constant is of order e, such that close to the critical 
dimension the RG expansion is under control. This flxed 
point controls the large-scale behavior of equation (||) for 
d < 5 and is associated with the non-trivial exponents 

The scaling exponent x of the interface height is zero 
for d = 2, while it is negative, x = —1/5, for d = 3. 
From equation (|^), with xo and zq replaced by the new 
exponents x and z, we see that in the two-dimensional 
case the highly elongated domains are only marginally 
stable against thermal fluctuations, since L|j- = 0(1), 

^x/z — 0(1), and L± = 0(1) (up to logarithms in 
each case). A thermally-induced stretching-and-breaking 
mechanism cannot, therefore, be excluded in this case. 
On the other hand, for d = 3 a negative value of x 
means that thermal fluctuations of the interfaces satu- 
rate at late times. Thus, in the three-dimensional case 
thermal roughening is depressed by the shear and the 
flow-induced layered structure of the domains is stable 
against thermal fluctuations. 

In this work we have shown that interface fluctuations 
under shear, in a system with non-conserved order pa- 



rameter, can be described by a stochastic differential 
equation, with an anisotropic Burgers non-linearity, and 
that the critical exponents are known exactly. In this 
way it was possible to assess the stability of domains in 
the late time regime of a system subject to a uniform 
shear flow. 

AB thanks R.K.P. Zia for a useful discussion. This 
work was supported by EPSRC grant GR/L97698 (AC 
and AJB), and by Fundagao para a Ciencia e a Tecnolo- 
gia grant BD/21760/99 (RDMT). 



[1] T. Hashimoto, K. Matsuzaka, E. Moses, and A. Onuki, 
Phys. Rev. Lett. 74, 126 (1995); J. Lauger, C. Laubner, 
and W. Gronski, Phys. Rev. Lett. 75, 3576 (1995). 

[2] A. J. Wagner and J. M. Yeomans, Phys. Rev. E 59, 4366 
(1999); 

[3] F. Corberi, G. Gonnella, and A. Lamura, Phys. Rev. 

Lett. 81, 3852 (1998); N. P. Rapapa and A. J. Bray, 

Phys. Rev. Lett. 83, 3856 (1999). 
[4] A. J. Bray and A. Cavagna, J. Phys. A 33, L305 (2000); 

A. Cavagna, A. J. Bray, and R. D. M. Travasso, Phys. 

Rev. E 62, 4702 (2000). 
[5] T. Ohta, H. Nozaki, and M. Doi, Phys. Lett. A 145, 304 

(1990); J. Chem. Phys. 93, 2664 (1991). 
[6] F. Corberi, G. Gonnella, and A. Lamura, Phys. Rev. 

Lett. 83, 4057 (1999). 
[7] S. F. Edwards and D. R. Wilkinson, Proc. R. Soc. Lon- 
don, Ser. A 381, 17 (1982). 
[8] J. M. Burgers, The nonlinear Diffusion Equation (Reidel, 

Boston, 1974). 

[9] M. Kardar, G. Parisi, and Y.-C. Zhang, Phys. Rev. Lett. 
56, 889 (1986). 

[10] E. Medina, T. Hwa, M. Kardar, and Y.-C. Zhang, Phys. 

Rev. A 39, 3053 (1989). 
[11] T. Hwa and M. Kardar, Phys. Rev. A 45, 7002 (1992). 

Note that the dimension d in this paper is our d' . 
[12] V. Becker and H. K. Janssen, Phys. Rev. E 50, 1114 

(1994). 

[13] D. Forster, D. R. Nelson, and M. J. Stephen, Phys. Rev. 
A 16, 732 (1977). 



3 



